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ABSTRACT 

The GOES M2-class solar flare, SOL2010-06-12T00:57, was modest in many 
respects yet exhibited remarkable acceleration of energetic particles. The flare 
produced an ~50 s impulsive burst of hard X- and 7-ray emission up to at least 
400 MeV observed by the Fermi GBM and LAT experiments. The remarkably 
similar hard X-ray and high-energy 7-ray time profiles suggest that most of the 
particles were accelerated to energies >300 MeV with a delay of ^10 s from mildly 
relativistic electrons, but some reached these energies in as little as s. The 
7-ray line fluence from this flare was about ten times higher than that typically 
observed from this modest GOES class of X-ray flare. There is no evidence 
for time-extended >100 MeV emission as has been found for other flares with 
high-energy 7 rays. 

Subject headings: Acceleration of particles — Sun: flares — Sun: particle emis- 
sion — Sun: X-rays, gamma rays 


1. Introduction 

The Sun is capable of accelerating electrons and ions to relativistic energies on time 
scales as short as a few r seconds during solar flares. This conclusion has been reached based 
on observations of the X rays, microwaves, 7 rays, and neutrons produced when the flare- 
accelerated particles interact in the solar atmosphere (Forrest & Chupp 1983; Kane et al. 
1986). The first reported observation of 7 rays with energies above 10 MeV was made with 
the Solar Maximum Mission ( SMM) spectrometer during the 1981 June 21 flare (Chupp et al. 
1982). Most of the 7-ray emission occurred within an ~70 s period and was followed minutes 
later by detection of high-energy neutrons at the spacecraft. Although it was clear from the 
neutron timing observations that protons were accelerated to energies in excess of 100 MeV, 
it was not possible from 7-ray spectroscopic studies to conclude that the protons reached the 
energies above ~300 MeV necessary to produce the characteristic spectrum from pion-decav 
radiation (Murphy et al. 1987). 
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Forrest et al. (1985, 1986) provided clear spectroscopic evidence for pion-decay emission 
during both the prompt and delayed emission phases of the 1982 June 3 flare. This emission 
in the second phase was confirmed by the presence of the 0.511-MeV annihilation line from 
the decay of positively charged pions that had a time profile similar to the high-energy 7-rays 
(Share et al. 1983). 

Since the era of those early measurements, improved spectrometers have detected the 
presence of pion-decay emission in several more flares and in two cases have observed 7-ray 
emission up to at least 1 GeV (Kanbach et al. 1993; Akimov et al. 1996; Vilmer et al. 2003; 
Kuznetsov et al. 2011). Ryan (2000) and Chupp & Ryan (2009) have reviewed these and 
other examples of high-energy 7-ray emission in flares. Most of these high-energy emissions 
have been observed over tens of minutes to hours leading to the designation of a class known 
as ‘Long Duration Gamma-Ray Flares’. Several processes have been suggested to explain 
such particle acceleration to high energies in the solar environment (Ellison &; Ramatv 1985; 
Petrosian et al. 1994; Park et al. 1997; Ryan 2000; Aschwanden 2004; Chupp & Ryan 2009). 

Typically such high-energy 7-ray emission has been associated with intense soft X-ray 
flares recorded by GOES' having peak powers exceeding 10" 4 W m“ 2 (X class). It is of interest 
to determine whether less energetic solar flares have the capability of accelerating electrons 
to energies >100 MeV and protons to energies >300 MeV. Out of 24 flares observed at > 10 
MeV energies by the 7-ray spectrometer on SMM only 3 had GOES classifications of M5 or 
less (Vestrand et al. 1999). Similarly, most flares emitting 7-ray lines have been associated 
with intense soft X-ray emission: out of 65 7-ray line flares observed by SMM only three have 
been associated with flares having a GOES classification of M5 or less. Before detection of 
this flare, the Ramaty High Energy Solar Spectroscopic Imager ( RHESSI) observed only one 
nuclear-line flare out of 20 below r the M5 class (Shih et al. 2009) and that was characterized 
as C9.7. The smallest GOES class flare for which detection of nuclear 7-rays has been 
claimed was the C7 flare observed by the COMPTEL Compton telescope on the Compton 
Gamma-Ray Observatory (CGRO) (Young et al. 2001). 

With the launch of the Fermi mission in 2008, it is now possible to make the high- 
sensitivitv measurements necessary to detect > 30 MeV 7 rays in the weakest flares. Fermi 
is comprised of two instruments: the Gamma-ray Burst Monitor (GBM) (Meegan et al. 
2009) sensitive from ~ 8 keV up to 40 MeV and covering the energy band of nuclear 7-ray 
line emission; and the Large Area Telescope (LAT) (Atwood et al. 2009) operating from 20 
MeV to more than 300 GeV and covering most of the pion-decay emission energy range. In 
its normal sky survey mode the LAT observes the Sun for only ~ 35 min every 3 hours. It 
was therefore fortunate to observe the M2 GOES-cl&ss flare on 2010 June 12 (SOL2010-06- 
12T00:57). This flare produced nuclear line emission about an order of magnitude higher 



than typical for this magnitude soft X-ray flare. Although the flare lasted only about 2 
min, it appears to have accelerated electrons and/or protons to energies >300 MeV based 
on detection of 7 rays with energies >100 MeV. One of the key features of this observation 
is that the high-energy emission is delayed by ~ 5 — 10 s relative to the 500 - 1000 keV 
bremsstrahlung. This suggests that the acceleration of hundreds of keV electrons and hun- 
dreds of MeV electrons and/or protons in the chromosphere took place within 10 s of each 
other. A weak solar energetic particle (SEP) event observed by GOES followed this flare. 

In the next two sections we discuss the capabilities of the GBM and LAT for detecting 
7-rays from solar flares. These are followed by §4 describing the observations including 
spectroscopic studies in the nuclear energy range with the GBM and in the high-energy 
domain with the LAT. As most high-energy 7-ray flares previously observed exhibited an 
extended emission phase, we then discuss a search for such emission in this flare in the next 
section. In the final section, §5, we summarize our conclusions and discuss the results. 


2. GBM Capability for Flare Observations 

The Gamma-ray Burst Monitor (GBM) was designed to observe gamma-ray bursts but 
has useful capabilities for other sources such as SGRs (soft 7-ray repeaters), pulsars, X-ray 
binaries using the Earth occultation technique, Terrestrial Gamma-ray Flashes (TGFs), and 
solar flares (Meegan et al. 2009). The GBM is comprised of twelve sodium iodide (Nal) 
detectors measuring the energy range from keV to 1 MeV, and two bismuth germanate 
(BGO) detectors covering the range from ~200 keV to 40 MeV. The detectors are arranged 
to collectively view the entire sky not blocked by the Earth. 

In this paper we use the BGO detector viewing the Sun to observe the nuclear line 
and continuum emission from solar flare SOL2010-06-12T00:57. Each BGO detector is a 
cylinder of length and diameter 12.7 cm, viewed by a photomultiplier tube at each end. Its 
effective area for detecting photons ranges from 160 cm 2 to 200 cm 2 , depending on energy 
and direction of incidence. Importantly for measuring solar flare nuclear lines, the high Z of 
bismuth, the high density of BGO, and the large volume of the GBM BGO detectors result 
in a high probability for absorbing the full incident photon energy: ^ 67% at 1 MeV, 50% 
at 3 MeV and 40% at 10 MeV. 

Several types of data are produced by the GBM. There are two temporally-binned types, 
Continuous Time (CTIME) and Continuous Spectroscopy (CSPEC), as well as Time-Tagged 
Events (TTE). For intense solar flares, TTE data are likely to be lost due to a bandwidth 
limit between the GBM and the spacecraft solid-state recorders. Therefore for solar-flare 
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spectral analysis the appropriate data type is CSPEC. CSPEC has 128 quasi-logarithmic 
energy bins and 4.096 s temporal resolution, which is improved to 1.024 s wdien the flight 
software triggers on a statistically significant rate increase. 

The detector performance was measured before launch using X-rays from the BESSY 
synchrotron radiation facility, radioactive isotopes including the 12 C line at 4.43 MeV (from 
an 241 Am/ 9 Be source), and 7 ~rays up to 17.6 MeV produced via (p, 7 ) reactions using a 
Van de Graaff accelerator. These data have been used to calibrate the gain, resolution, 
and response of the detectors (Bissaldi et al. 2009; Meegan et al. 2009). The automatic 
gain software measures and adjusts the gains of the BGO detectors using the 2.223 MeV 
line present in the background from capture of cosmic-ray produced neutrons on H in the 
spacecraft’s propellant. In order to improve the gain solution for analysis of the June 12 
flare data we measured the centroid of the solar neutron-capture line. This resulted in a 
1% gain change that is small relative to the ^ 8% resolution of the GBM at the 2.223 MeV 
line. The best-fit energy of the positron annihilation line (0.511 MeV) from the flare was 
found to be 0.530 ± 0.007 MeV, in disagreement with fits of the instrumental background line. 
However, both lines are superimposed on strong continua making it difficult to determine the 
lines 3 * * * 7 peak position. Further investigation of the BGO calibration below 1 MeV is planned, 
particularly since TGF spectra also show a shifted positron annihilation line (Briggs et al. 

2011). 

Spectral fitting is performed using the forward- folding technique with an assumed parametrized 
photon spectrum folded through a detector response matrix (DRM) to produce a counts 
spectrum; the parameters are adjusted to obtain the best fit to the observed counts spec- 
trum. The DRM is based on Geant4 simulations of the GBM detectors and mass model of 
the satellite and has been validated by comparison with the calibration data (Hoover et al. 

2008). 


3. LAT Capability for Flare Observations 

The Large Area Telescope (LAT) is a pair-conversion telescope designed to detect 

gamma rays from 20 MeV up to more than 300 GeV (Atwood et al. 2009). It is made up of a 
4x4 matrix of identical towers, each one comprised of a tracker with layers of Silicon Strip 

Detectors (SSD) alternating with foils of high-Z converter (tungsten), and a calorimeter with 
logs of Csl arranged in a Todoscopic 7 configuration so that the energy deposition is imaged 
in three dimensions. The array of towers is surrounded by an Anti-Coincidence Detector 

(ACD) made up of 89 plastic scintillator tiles with a 5 x 5 array on the top and 16 tiles 
on each of the 4 sides. To reduce the impact of the self- veto events caused by calorimeter 



back-splash, each ACD tile is typically <1000 cm 2 , depending on its position in the array. 
A 7-ray passes through the ACD with small probability of interaction and can convert into 
an electron-positron pair that is tracked in the SSD. The energy of a photon below about 
100 MeV can be estimated by multiple-Coulomb scattering of the electrons in the tungsten 
converters and SSD. Energies of higher-energy photons are measured by ‘total’ absorption 
or modeling the shower profile in the calorimeter. In the standard Fermi sky-survey mode 
the spacecraft rocks 50° north and south in celestial declination from the zenith so that each 
region of the sky is viewed for ^30-40 min every two orbits; therefore the LAT’s large aper- 
ture (2.4 sr) and effective area provides the capability to sensitively monitor solar activity 
with a duty cycle of 15-20%. 

The ACD is used to reject the large background of charged cosmic-radiation and secon- 
daries from the spacecraft and Earth’s atmosphere. The threshold for this veto is nominally 
set at 45% of the amplitude of a minimizing ionizing singly-charged particle traversing a tile 
in the ACD (i.e. set to ~800 keV). If this threshold is exceeded and if the ACD tile hit 
was adjacent to a tracker tower that caused the event to trigger, the event will be vetoed 
unless a sufficiently high energy is deposited in the calorimeter (>100 MeV in one or more 
crystals). This veto inhibition ensures that backsplash does not cause very high-energy 7 
rays to self-veto. Of the vetoed events, 2% are telemetered to the ground for diagnostic 
purposes. 

Solar flares can emit intense fluxes of tens of keV X rays. About 20% of the X rays at 
these energies can penetrate the thermal blanket and micrometeoroid shield and reach the 
ACD, depositing part of their energy in the illuminated tile. 

Several of these photons can arrive within the 0.4 fis anti-coincidence veto shaping time 
(pulse pile-up) to yield a high total-energy loss. It is possible in the largest of flares that this 
energy loss could exceed the 800 keV veto threshold and information about a valid photon 
event would not be transmitted to the ground. None of the solar flares detected to date have 
been sufficiently intense to cause this data loss to occur. 

Pulse pile-up of hard X rays from some flares has, however, affected the ground analysis 
of LAT data, in the classification of events as 7 rays or background. ACD tile hits are 
registered (i.e. included in an event) whenever there is >100 keV energy deposition integrated 
over the ^3 (i s peaking time of the front-end amplifier (as opposed to the 800 keV threshold 
for an ACD veto). Pulse pile-up during even modest flares generates such tile hits. In order 
to achieve the highest sensitivity for studying celestial sources, the current LAT ground- 
processing software rejects 7-ray events with a high ratio between the number of ACD tile 
hits (i.e,, the number of tiles with energy deposition >100 keV) and the energy measured by 
the calorimeter. Thus, 7 rays arriving during relatively intense flares have a high likelihood 
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of being rejected from the standard LAT data products. This is the case for the 2010 June 
12 flare as we will discuss in §4.2 below. 

As a consequence of this feature, the LAT Low-Energy (LLE) technique (Pelassa et al. 
2010) was adapted for analyzing data during this flare. LLE event selection uses less dis- 
criminating criteria than the standard ground- processing and is not affected by ACD tile hits 
>100 keV. The primary requirement is that the candidate 7-rav event have at least one track 
and a reconstructed energy larger than 30 MeV. Due to these less discriminating criteria the 
off-axis angle for LLE events can be as large as ^80° compared with M)8° in the standard 
LAT data products. In addition, only 7-rays whose estimated arrival directions were within 
20° of the Sun are included for analysis. The 20° angular restriction was chosen based on 
Monte Carlo simulations of the LAT point-spread function at these energies. A time series 
of events is then constructed from which background intervals on either side of the flare are 
defined and a linear or quadratic interpolation is used to estimate the background during the 
flare. A Detector Response Matrix (DRM) for the solar location during the flare is created 
using a custom Monte Carlo simulation. By passing candidate photon models through the 
DRM, we then fit the background-subtracted data using a version of rmfit 3.4 1 , customized 
for the specific solar flare, and the OSPEX 2 analysis packages. 

Using a Monte Carlo analysis, we estimate that the energy resolution for a source at ~75° 
off axis in our LLE analyses is about 40%. We also estimate ^15% and ~10% systematic 
uncertainties in the LLE fluxes >30 and 100 MeV, respectively, based on studies of the Vela 
pulsar. The uncertainty >100 MeV is consistent with that found using the standard Fermi 
analysis tools applied to the pulsar (Abdo et al. 2010). 


4. June 12 Flare Observations 

The GOES M2.0 class X-ray flare commenced with some low-level activity on 2010 
June 12 at 00:30 UT. Although both RHESSI and Fermi observed the flare, here we only 
discuss the Fermi observations because RHESSI was not pointed at the Sun during the 
flare. The flare occurred with Fermi in sunlight and during a relatively low-background 
portion of its orbit. As seen in Figure la the 10 — 26 keV emission recorded by the GBM 
Nal detectors commenced around 00:55 UT and rose precipitously about 40 s later; for 
comparison we also plot the GOES 0.5 - 4 A profile and note that this emission is dominated 


1 R.S. Mallozzi, R.D. Preece, & M.S, Briggs, “RMFIT, A Lightcurve and Spectral Analysis Tool”, Robert 
D. Preece, University of Alabama in Huntsville, (2008): http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 

2 SolarSoffc: http://www.lmsal.com/ solarsoft / 
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by 3 keV thermal photons as is reflected in its slower rise and extended tail. The GBM burst 
algorithm triggered on the flare at 00:55:05.64 UT and put the instrument in a high-data 
rate mode for the next 10 min. The 100 — 300 keV time profile observed by the GBM’s solar 
facing Nal detector is also plotted in Figure la. It is clear that the emission peaks more 
sharply and ends sooner at higher X-ray energies. Most of the emission observed above 100 
keV occurred within an ~50 s interval. The event as viewed by GOES 1-8 A channel 
ended about 01:30 UT. The flare originated from active region (AR) 11081 at approximately 
N23°W43°. White light emission observed by the Helioseismic and Magnetic Imager (HMI) 
on the Solar Dynamics Observatory (SDO) (Martinez Oliveros et al. 2011) in a single 45 s 
exposure, consistent in time with the hard X-ray emission, reveals two compact footpoints 
about 10 4 km apart. 

There was no evidence for significant dead time and/or pulse pile-up effects in the 
GBM BGO detector facing in the solar direction. Photons with energies up to MeV 
were detected by the GBM during the 50 s peak. At the time of the flare the spacecraft was 
rocking 50° to the south so that the Sun w r as 76° off axis, close to the edge of the field-of-view 
(FOV) for LLE studies, and the Earth’s horizon was entering the FOV. The accompanying 
hard X-ray emission from the flare was detected in the LAT’s ACD. In Figure lb we plot the 
average number of ACD tile hits as a function of time. As discussed in the previous section, 
pulse pile-up from tens of keV hard X-rays exceeded the 100 keV ACD hit threshold. This is 
reflected in the broad peak with a maximum near 00:57 UT that has a shape similar to the 
11-26 keV emission observed by the GBM Nal detector. The impulsive peak in the ACD rate 
is also similar to that observed between 100 and 300 keV by the GBM Nal detector. There 
is no evidence for an increase in the number of ACD vetoed events in orbit during the flare, 
indicating that pulse-pileup from hard X-rays did not exceed the 800 keV veto threshold. 
Thus, the overall valid event rate transmitted to the ground for processing was not affected 
by the ACD response. However, as shown by the red curve in Figure lc there is no evidence 
for the flare in the well-screened standard LAT data products (shown in the figure are the 
events belonging to the “transient” event class, Atwood et al. 2009). If anything, we see a 
deficit of events in the standard analysis light curve which is a consequence of the high ACD 
hit rate >100 keV, shown in Figure lb, that caused a significant increase in the number of 
events that failed the standard quality cut. 

It is important to establish convincingly that any high-energy emission observed by the 
LAT originated at the flare site and was not due to artifacts from the high rates encountered 
by the instrument. The most compelling evidence for the solar flare origin is the map of 
events relative to the position of the Sun. Plotted in Figure 2a is the distribution of LLE- 
selected events with energies >30 MeV accumulated 30 s before and after the flare. The 
distribution is affected by the instrument field of view and by the removal of events near the 




Time (UT ? 2010-06-12) 

Fig. 1.— Time histories related to the 2010 June 12 solar flare, a) GOES 0.5 -4 A rates, 
and GBM Nal 11-26 keV and 100 - 300 keV relative rates; b) LAT ACD hit rate >100 keV 
containing contributions from background, >100 keV solar flare X rays (impulsive peak) and 
pulse pile up from 10's of keV solar X rays following the Nal 11-26 keV profile in la); and 
c) LLE and LAT Transient Class event rates. 
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Earth’s horizon. We reduce the contamination from the bright Earth limb selecting only the 
events with a reconstructed zenith angle less than 100°, but, due to the large point spread 
function at low energy, some residual contamination near the horizon may be due to this 
limb brightening. 

In Figure 2b we plot the distribution accumulated from 00:55:40 to 00:56:30 UT dur- 
ing the flare. We note that the excess near the Sun’s position is biased due to its loca- 
tion near the edge of the LAT’s field of view. This is known as the “fisheye” effect, i.e., 
the tendency of events to be reconstructed with directions closer to the z axis than they 
should (Thompson et al. 1993). This bias is relevant only for off-axis events, and it is par- 
ticularly evident at low energies. The position of the centroid of the counts is shown by the 
filled star (at R.A.=84°.39, Dec. =18°. 79, J2000). It is important to note, that in the routine 
LAT analyses the contribution to the overall exposure and photon counts from far off-axis 
events is negligible. In addition to this, owing to the scanning mode of observations with 
the LAT, persistent sources will have the great majority of their exposure at much smaller 
inclination angles. 

In Figure 2c we plot the background-subtracted distribution observed during the flare. 
There appears to be significant 7-ray emission consistent with the reconstructed position 
of the Sun up to energies of ~400 MeV. We have confirmed the observed 6° shift in the 
position of the centroid from the solar position using a Monte Carlo study. Using the shift 
and identifying a < 20° acceptance angle (dashed curve in the figure) for LLE events >30 
MeV maximizes the signal/background. Events that meet these criteria make up the LLE 
rate (black curve) plotted in Figure 1c. The overall rate decreased, especially after 01:00 
UT because the defined aperture approached the exclusion region near the Earth’s limb and 
many events were consequently rejected. The >30 MeV LLE time profile during the flare is 
similar to the 100-300 keV Nal time history. 


4.1. GBM Spectroscopic Studies 

In order to obtain background-subtracted spectra we used GBM/BGO spectra accu- 
mulated just before the flare and after 4-min following the flare for background in order to 
avoid times when there was significant 15-50 keV hard X-ray emission. The background- 
subtracted GBM counts spectrum accumulated over a 50-see period (00:55:40-00:56:30 UT) 
during the impulsive phase of the flare is shown in Figure 3, Although the flare was relatively 
weak and the BGO detector has only moderate energy resolution, line features are clearly 
evident in the spectrum. 
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R.A. 


Fig 2. — Angular distribution of >30 MeV 7 rays relative to the position of the Sun detected 
using the LLE analysis. The open star shows the location of the Sun wdiich is close to the is 
close to the edge of the LAT field of view. The filled star shows the shifted position due to 
the fisheve effect discussed in the text. The dotted curve depicts the 20° region containing 
the events used in the analysis, a) Average of background distribution taken 30 s before and 
after the flare, b) Angular distribution of events observed between 00:55:40 and 00:56:30 
UT. c) Difference between angular distribution observed during the flare and the average 
background distribution. The distributions have been smoothed (with a Gaussian kernel) to 
reduce statistical fluctuations. 
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Fig. 3. — Background-subtracted BGO count spectrum accumulated between 00:55:40 and 
00:56:30 UT. The spectrum has been fit by a simple power law (blue histogram) and a 
flatter power-law with exponential cut off (light green histogram) to model bremsstrahlung 
by electrons (see text), and with a nuclear de-excitation component (purple) plus 0.511 and 
2.223 MeV lines (dashed and dark green), and a pion-decay (orange) component. 
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The spectrum has been fit with a photon model consisting of several components. We 
find that the electron-bremsstrahlung component has the shape of power law at low energies 
that hardens above several hundred keV and then rolls over in the MeV range. In this paper 
we fit the bremsstrahlung spectrum with the sum of a low-energy power law and a fiat power 
law with exponential cutoff in the MeV range (blue and light green curves, respectively, in 
Figure 3). Such a complex shape has been observed before in the spectra of several flares 
detected by SMM and RHESSI (Share &; Murphy, in preparation). McTiernan & Petrosian 
(1990) have found that the magnitude of the hardening above a few hundred keV observed in 
some flares is larger than that expected for an electron spectrum following a single power law. 
Rieger Sz Marschhauser (1991) describe flattening in the MeV range followed by roll overs 
above several tens of MeV in the spectra of some flares. Such features cannot be explained 
by transport effects alone (Petrosian et al. 1994) and must be produced by the acceleration 
mechanism. Park et al. (1997) show that these features can be explained by models based on 
stochastic acceleration by turbulence once loss mechanisms are properly included. Whether 
such models can explain the 2010 June 12 bremsstrahlung spectrum requires more study. 

The nuclear de-excitation lines and continua are represented by a template based on a 
detailed study of nuclear gamma-ray production from accelerated-particle interactions with 
elements found in the solar atmosphere (Murphy et al. 2009). Such templates depend on 
the assumed ambient composition and accelerated-particle composition and spectrum. The 
BGO spectral data for this flare are inadequate to distinguish among templates derived for 
different ambient abundances, particle spectra and angular distributions. For this reason 
we used a arbitrary templates based on earlier studies (Share &; Murphy, in preparation). 
The accelerated particles were assumed to interact in a thick target with a coronal compo- 
sition (Reames 1995), but with 4 He/H = 0.1. We do not take into account transport effects 
but instead assume that the accelerated particles have a power-law differential spectrum 
(dN/dE oc with f3 = —4), coronal elemental abundance (but with an accelerated a/p 
ratio of 0.2), and an angular distribution that is isotropic in the downward hemisphere. As 
we mentioned above our results are not sensitive to the assumed compositions. The pho- 
ton model also includes Gaussians representing the 0.511 and 2.223 MeV positron-electron 
annihilation and neutron-capture lines, respectively, and a pion-decay spectral component. 

We present the best- fitting spectral parameters in Table 1 , along with estimates of their 
la uncertainties. We list fluences obtained by integrating over the 50-s time period, except 
for the 0.511 and 2.223 MeV lines. The 0.511 MeV line originates from radioactive decays 
(Kozlovsky et al. 1987, 2004) with half lives that can extend up to hours, in addition to a 
prompt component from positron annihilation following decay of positively charged pions. 
The neutron-capture line is also delayed because of the time required for the neutrons to slow 
down in the solar atmosphere and photosphere and be captured on H, forming deuterium 
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with the release of 2.223 MeV of binding energy from the mass excess. For these two lines 
we have fit the background-subtracted spectrum integrated over a total of 250 s. 

The low-energy power law r (PL1) is well defined by the fit. The higher-energy power law 
(PL2) with exponential rollover is not as well defined because it competes with the nuclear 
de-excitation line spectrum. However our studies indicate that this component is required 
to provide an acceptable fit to the June 12 flare spectrum. With this component present, 
the fit is acceptable at the 10% confidence level based on the y 2 statistic; without it, the fit 
is rejected (0.001 % confidence level). There is no evidence for pion-decay emission in the 
BGO spectrum. As can be seen in Table 1, assuming an hadronic origin, the LAT-measured 
pion-decay fluence that we discuss in the next section is about a factor of ten lower than 
detectable by the GBM BGO detector. 


4.2. Combined LAT and GBM Spectroscopic Studies 

We have obtained a background-corrected LAT count spectrum >30 MeV accumulated 
during the 50-s period 00:55:40-00:56:30 UT using the LLE data plotted in Figure lc. This 
spectrum revealed flare emission up to an energy of ~400 MeV. The fundamental question 
is: what is the origin of this emission? The nuclear line emission observed with the GBM 
implies the presence of accelerated ions up to at least 50 MeV nucleon -1 . It is possible that 
the flare- accelerated proton spectrum extended up to the ~300 MeV threshold for pion pro- 
duction. Alternatively, it is possible that the LAT emission is from electron bremsstrahlung, 
either from an extension to high energies of the electron spectrum producing the X-ray 
bremsstrahlung observed in the GBM or from an additional hard electron component. One 
possible way to resolve this ambiguity is to jointly fit the GBM and LAT spectra assuming 
different origins for the LAT emission. 

In Figure 4 w r e plot the background-subtracted photon spectrum from 0.3 to 400 MeV 
including both the GBM and LAT data. We made two fits to the joint data sets, one 
assuming that the observed LAT emission was from pion-decay radiation (top panel) and 
the other assuming that it was from a hard power-law spectrum of electron bremsstrahlung 
(bottom panel). Based on the statistical quality of the fits to the LAT spectrum we cannot 
distinguish between the two emission models. In addition we cannot constrain the origin of 
the emission for this event by extrapolating the models into the GBM energy range; however 
we note that for a stronger flare we might be able to rule out a power-law model. The choice 
of model also has only a small effect on the parameters listed in Table 1 derived from fits to 
the 0.3 to 40 MeV GBM data. Also plotted in the figure are the extensions into the LAT 
energy range of the power-law and cutoff power-law components derived from the fits to the 
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Fig. 4— Combined GBM/LAT photon spectrum accumulated between 00:55:40 and 
00:56:30 showing the best total fit using the same components as in Figure 3 plus an addi- 
tional component for the LAT emission. The upper panel shows a pion-decav fit to the LAT 
spectrum; alternatively the lower panel shows a power-law fit, presumedly representing a 
third electron bremsstrahlung component. Note that because this is a photon representation 
the lines are plotted at their intrinsic resolution and appear to be more significant than they 
really are. 
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GBM data. The intensities of these components fall at least an order of magnitude below 
the LAT measurements and therefore do not make a significant contribution to the solar 
emission observed by the LAT. 

Even though we cannot statistically distinguish between a pion-decay or electron-bremsstrahlung 
origin for the observed LAT emission, we can obtain the best- fitting parameters for these 
components. If the LAT emission is from electron bremsstrahlung, we have shown that it 
cannot be a simple extension of the low-energy bremsstrahlung components that we deter- 
mined from fits to the GBM data; it must be from a distinct population of electrons extending 
to energies of several hundred MeV. We have determined the power-law parameters (PL3) 
of the fitted bremsstrahlung spectrum and list them in Table 1. However, this third high 
energy electron component cannot be produced by the acceleration models mentioned §4.1 
above, which produce spectra that steepen beyond tens of MeV due to synchrotron energy 
losses that increase with energy (see Park et al. 1997), and must have a quite different origin. 
Consequently we believe that this is a less likely scenario than the hadronic model. 

Assuming that the LAT emission is from hadronic interactions, we have fit the LAT 
spectrum with calculated pion-decay spectra produced by accelerated ions having differential 
power-law indices from —2.5 to —7.5. With 67% confidence (based on x 2 ) we conclude that 
the spectrum of accelerated ions responsible for the pion-decay emission must be steeper 
than a power-law with index —4.5. We note that there is no change in the quality of the 
fits for indices steeper than —5 due to limited statistics >400 MeV. We list the fluence of 
pion-decay photons >200 keV derived from these fits in Table 1. This fluence is a factor of 
ten below the limit derived from fitting the GBM BGO data up to 40 MeV, illustrating the 
significantly greater sensitivity of the LAT instrument. 

We can use the results of our GBM and LAT spectral analyses to obtain information 
on ions accelerated in the impulsive phase of the June 12 flare. Murphy et al. (1997) have 
described how r parameters derived from integrated spectroscopic fits and temporal studies 
can be used to obtain this information. We first use the nuclear de-excitation line, neutron- 
capture line, and pion-decay fluences listed in Table 1 to estimate the overall shape of the 
accelerated ion spectrum. These three emissions are produced by accelerated ions within 
distinct energy ranges: ~5-20 MeV for the de-excitation lines, ~ 10-50 MeV for the neutron 
capture line, and >300 MeV for the pion-decay emission. Ratios of these emissions therefore 
determine the relative numbers of accelerated ions in the associated energy ranges. We then 
obtain spectral indices across these energy ranges by comparing measured ratios with ratios 
from theoretical calculations (Murphy et al. 1987, 2005, 2007) based on updated nuclear 
cross sections. 

If we assume that the LAT emission >30 MeV was entirely due to pion-decay emission, 



-19- 


then we estimate that the flare- accelerated ion spectrum was consistent with a series of power 
law T s, softening with energy, with indices of 3.2 between ~ 5 — 50 MeV, ^—4.3 between 
^50-300 MeV, and softer than ~™4.5 above 300 MeV. However, these calculations assume 
that these individual power laws continued without break to high energies, which cannot be 
the case wdiere the spectrum softens with increasing energy. We will describe a more refined 
representation of the accelerated proton spectrum in a forthcoming publication. 


4.3. Combined GBM and LAT Timing Studies 

The combined GBM and LAT observations also provide us with the opportunity to 
study acceleration and transport in this impulsive flare. In Figure 5 we plot 5-s resolution 
time profiles of the fitted bremsstrahlung and nuclear de-excitation line fluxes from the GBM 
and the >100 MeV flux observed by the LAT. For purposes of comparison we have plotted 
the bremsstrahlung profile over the nuclear and LAT >100 MeV histories in b) and c), 
respectively. We have also plotted the nuclear profile over the LAT >100 MeV history. The 
early peaking of the bremsstrahlung suggests that the higher-energy emissions were delayed 
by a few seconds and a cross-correlation study indicates that the overall lag between the 
LAT >100 MeV flux and the GBM 300 keV bremsstrahlung flux is ~3 s. 

This delay warranted a higher-time resolution comparison of 100 - 500 keV GBM and 
LAT >30 MeV counting rates rather than fitted fluxes. The 100 - 500 keV band is dominated 
by electron bremsstrahlung as can be seen in Figure 3. In Figure 6 we plot GBM/BGO 
100-500 keV rates at 320 ms resolution with 3-s LAT/LLE >30 MeV rates over plotted. 
The hard X-ray profile reveals the presence of a clearly separated initial peak along with 
other structures. The onset of the > 30 MeV emission appears to be ~3 s following the 
bremsstrahlung and rises to a peak about 10 s after the 100 -500 keV peak. The LLE 
profile appears to reflect the double-peaked bremsstrahlung profile w r ith a delay of about 10 
seconds. From a cross correlation analysis of the two profiles plotted in Figure 6 w T e find 
that the >30 MeV emission lags the bremsstrahlung by 6 ± 3 s. There are two fundamental 
implications of the time profiles in Figure 6: 1) protons and/or electrons began reaching 
energies above 100 MeV within a few seconds of the time it took to accelerate electrons to 
energies of hundreds of keV; and 2) the overall acceleration time scale of the >100 MeV 
particles is similar to that observed in hundreds of keV electrons, but delayed by about 10 
seconds. 



Table 1. Best-fitting spectral parameters 


Parameter Value 

2.85 ± 0.1 7 cm ' 2 keV ” 1 


PL1 fluence at 300 keV 
PL 1 index 

PL2 fluence at 300 keV 
PL2 index 

PL2 Exponential Energy 
0.511 MeV line fluence a 
2.223 MeV line fluence a 
Nuclear line fluence 
Pion-decay fluence (GBM) > 200 keV 
Pion-decay fluence (LAT b ) > 200 keV 
Pion-decay fluence (LAT) >100 MeV 
PL3 fluence at 30 MeV 
PL3 index 


3.31 ± 0.09 

0.08 ± 0.02 7 cm ' 2 keV ” 1 
£ 1-2 

2400 ± 800 keV 

11.3 dfc 2.5 7 cm ” 2 

21.3 ± 2.0 7 cm ' 2 

23.5 d= 2.5 7 cm ” 2 

1.5 db 2.5 7 cm ” 2 
0.62 db 0.07 7 cm ” 2 
0.13 ± 0.015 7 cm ” 2 

(9.2 db 2.0) xl 0~ 6 7 cm ” 2 keV ” 1 
1.9 ± 0.2 


a lntegrated from 00:55:40 to 00:59:50 UT 

b Computed by extrapolating to low energies the model that best fits LAT 
data 



Time, UT 

Fig. 5. — Time profiles of the a) bremsstrahlung differential flux at 300 keV (GBM), b) total 
nuclear de-excitation line flux (GBM), and c) integral pion-decay flux >100 MeV (LAT). 
Solid curves over-plotted in b) and c) are the arbitrarily normalized bremsstrahlung flux 
plotted in a); dashed curve over-plotted in c) is arbitrarily normalized nuclear line flux from 
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4.4. Search for >100 MeV Emission Following the Flare 

The Fermi LAT detects quiescent emission from the Sun on a near daily basis (Abdo et al. 
2011). This emission comes from cosmic-ray proton interactions in the solar atmosphere 
and photosphere, and from Compton scattering of cosmic-ray electrons on solar black- 
body photons. The LAT is therefore a sensitive monitor of temporally extended solar-flare 
emission such as detected by CGRO EGRET experiment following the 1991 June 11 flare 
(Kanbach et al. 1993; Rank et al. 2001). We therefore studied the emission within 15° of 
the Sun in the hours preceding and following the flare. The standard LAT data products 
were used in the analysis, which modeled the region around the Sun including all sources in 
the Second Fermi LAT Source Catalog (The Fermi-LAT Collaboration 2011), isotropic and 
Galactic diffuse emissions, and spatially extended Compton-scattered solar photons discussed 
above. 

In Figure 7 we plot 95% confidence limits on the >100 MeV flux from the solar disk in 
^30 min exposures every two orbits from 6 hr before the flare to 22 hr after it. We note 
that the Sun was outside the FOV for LAT standard-product analysis during the flare, but 
its 7-ray emission could be studied in the orbits just before and after the flare. During this 
time period the Moon passed within 10° of the Sun and its flux of ~ 1 x 10“ 6 7 cm” 2 s” 1 
made a significant contribution to the measured solar fluxes because we did not include it 
in the model of the source region. There is no evidence for an increase in the solar emission 
following the time of the flare, which is denoted by the peak flare flux >100 MeV. This peak 
flux is about 1000 times higher than the plotted upper limits. We wish to compare this to 
the time extended >100 MeV emission observed by EGRET following the 1991 June 11 flare 
(Kanbach et al. 1993; Rank et al. 2001). But EGRET was saturated at the peak of the flare 
and therefore its time history cannot be normalized to the >100 MeV LAT peak flux of the 
June 12 flare. However, after the peak of the flare the ratio of the >100 MeV EGRET flux 
to the 2.223 MeV neutron-capture line measured by COMPTEL, which operated normally 
throughout the June 11 flare, was relatively constant. We therefore used the COMPTEL 
time history as a proxy for the EGRET >100 MeV photon fluxes. The time profile of the 
COMPTEL 2.223 MeV neutron-capture line flux (Share et al. 1997), normalized to the peak 
>100 MeV LAT June 12 flux, is shown in the Figure 7. We see that during the first 30 
minute exposure following the flare, the LAT upper limit is a factor of ~ 20 below the value 
expected if there had been comparable extended emission similar to that found on 1991 June 
11 . 
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Fig. 6. — Comparison of high-time resolution profiles of the 100 - 500 keV emission ob- 
served in the GBM BGO detector and of >30 MeV LLE data. A cross-correlation analysis 
indicates that the high-energy 7-ray emission had an overall lag of 6 ± 3 s relative to the 
bremsstrahlung. 
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Fig. 7. — 95% confidence upper limits of >100 MeV solar 7- rays measured by LAT within 7 
hours preceding and 22 hours following the 2010 June 12 flare. Dashed line is the 1991 June 
11 time history of the 2.223 MeV neutron-capture line observed by COMPTEL experiment 
and diamonds are EGRET >50 MeV data both normalized to the peak >100 MeV flux 
measured by the LAT for the June 12 flare. The 2010 June 12 flare does not show extended 
emission like the 1991 June 11 flare observed with EGRET. 
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5. Summary and Discussion 

The 2010 June 12 flare was the first in Cycle 24 to be observed to emit nuclear 7 
rays. It was also the first flare detected by the Fermi LAT at energies above 30 MeV. The 
hard X-ray and nuclear line radiation was observed both by the Fermi GBM and RHESSI 
spectrometers. In this paper we only analyzed GBM data because RHESSI was offset from 
the Sun to study the Crab Nebula during the time period of the flare; this affected knowledge 
of the instrument response. 

The fact that the flare emitted detectable 7 rays at all is surprising because its peak soft 
X-ray emission only reached a GOES M2 level. However, Shih et al. (2009) and others have 
shown that 7-ray line fiuences are only weakly correlated with GOES soft X-ray emission 
but are strongly correlated with electron bremsstrahlung fiuences >300 keV. This is true for 
the June 12 flare as we find that the measured bremsstrahlung and 2.223 MeV fiuences are 
consistent with the established correlation. 

The flare originated from a compact region and its hard X-ray emission only lasted 50 
s. Figure 6 reveals striking information about the processes that accelerate protons and/or 
electrons to energies of hundreds of MeV. We find that although some of the particles reach 
energies >100 MeV within about 3 s, the bulk of these particles reach such high energies 
following a delay of about 10 s. This is revealed in the delayed double-peaked time structure 
>30 MeV that is similar to what is observed in hard X-rays. 

In Table 1 we list the best-fitting parameters from our fits to the GBM and LAT spectra. 
These include the amplitudes (at 300 keV) and indices of two power-law continua observed by 
the GBM between 300 keV and 8 MeV. The first is an extension of the hard X-ray spectrum 
observed by the GBM Nal detectors. The second appears to be a hard power law with an 
exponential cutoff energy near 2.5 MeV. Although the GBM only has moderate spectral 
resolution, it was able to measure the fiuences of the 0.511 MeV annihilation and 2.223 MeV 
neutron capture lines, and the total nuclear de-excitation emission. There was no evidence in 
the GBM data for flare emission above about 8 MeV. In contrast the LAT detected significant 
continuum emission from ^30 to 400 MeV about an order of magnitude below the GBM 
upper limits. This reflects the excellent sensitivity of the LAT for observing solar flares. 
This radiation could be either pion-decay emission or primary electron bremsstrahlung. In 
Table 1 we list the measured pion-decay fiuences >200 keV and >100 MeV. We also list the 
parameters of a high-energy power-law bremsstrahlung component that fits the LAT data 
equally well. Theoretical arguments made in §4.2 weigh against an electron origin for the 
emission observed by the LAT. Under the assumption of an hadronic origin for the LAT 
emission and using 7-ray line measurements by the GBM, we have estimated the shape of 
the accelerated-ion spectrum that could have produced the combined spectrum. 
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We have also set significant constraints on any time-extended >100 MeV emission. Our 
limit on the >100 MeV flux of photons is about an order of magnitude below what would 
have been expected if the decay followed that observed in the well-studied 1991 June 11 
flare (Kanbach et al. 1993; Rank et al. 2001). Thus, there is no evidence for precipitation 
of trapped flare particles, particles accelerated in magnetic loops after the impulsive phase, 
particles accelerated in CME-associated reconnection sheets (Ryan 2000), or particles sharing 
the same origin as the Solar Energetic Particles (SEPs) observed in space (Ramaty et al. 
1987; Oliver et al 1993). We estimate from white light measurements of the 1991 June 11 
flare (Sakurai et al. 1992) that the foot-point separation was ~ 2.5 x 10 4 km, about 2.5 
times larger than the June 12 flare. It is possible that longer coronal loops are necessary for 
time-extended acceleration and/or trapping of protons. 

This paper primarily addresses only the 7-ray observations of the 2010 June 12 flare. 
Other studies are currently in progress involving hard X-ray observations of the flare by the 
Fermi GBM Nal detectors and comparisons of the characteristics of the ion and electron 
populations at the Sun and in space and at Earth. 
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